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I ' ABSTRACT 

We present the results of an optical spectroscopic study of the massive Wolf-Rayet 
binary HD192641=WR137. These 1986-2000 data cover the dust-formation maximum 
in 1997. Combining all available measurements of radial velocities, we derive, for the 
first time, a spectroscopic orbit with period 4766 ± 66 days (13.05 ± 0.18 years). 
The resulting masses, adopting i=67°, are Mo=20 ± 2M Q for the O component and 
Mwr=4.4 ± 1.5M© for the WR component. These appear, respectively, ~ normal 
and on the low side for the given spectral types. Analysis of the intense multi-site 
spectroscopic monitoring in 1999 shows that the CIII A5696 and CIV A5802/12 lines 
have the highest intrinsic variability levels. The periodogram analysis yields a small- 
amplitude modulation in the absorption troughs of the CIV A5802/12 and Hel A5876 
lines with a period of 0.83 days, which could be related either to pulsations or large- 
scale rotating structures as seen in the WN4 star EZ Canis Majoris (WR6). Wavelet 
analysis of the strong emission lines of CIII A5696 and CIV A5802/12 enabled us to 
isolate and follow for several hours small structures (emission subpeaks) associated 
with density enhancements within the wind of the Wolf-Rayet star. Cross-correlating 
the variability patterns seen in different lines, we find a weak but significant correlation 
between the varability in emission lines with different ionization potential, i.e. in lines 
formed at different distances from the WR stellar core. Adopting a (3 wind-velocity 
law, from the motion of individual subpeaks we find [3 ~ 5, which is significantly larger 
than the canonical value f3 ~ 1 found in O-star winds. 

Key words: stars: Wolf-Rayet - stars: wind, outflows - stars: oscillations - binaries: 
spectroscopic 
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1 INTRODUCTION 

Population I Wolf-Rayet stars are highly evolved He-burning 
descendents of massive OB-stars, which are assumed to 
have undergone rapid luminous blue variable (LBV) mass- 
loss e pisodes before reaching the WR stage jLanger et alJ 
Il994l) . They exhibit intense b road emissio n lines pro- 
duced by highly-ionized atoms llConti 1 1 2000|) , which form 
a hot , fast and dense stellar wind dHamann fc K oestcrkc 
1998) . Furthermore, rec e nt observations iMoffat et alJll988t 
iLepine fc Moffat I ll99St ILepine et all l200d) revealed that 
the winds are highly structured, although the "Standard 
Model" of a WR-star wind was based on the simplify- 
ing assumptions of spherical symmetry, homogeneity, time 
independence and a monoton ic velocity law for the wind 
Jffillier Ill99~0: lHamann Ill99ll) . The highly structured na- 
ture of the winds may be related to numerous small- 
scale density enhancements embedded in the outmoving 
WR wind. Once included in the framework of the "Stan- 
dard Model", the micro-structuring (clumping) allowed one 
to: (i) produce much better fits of the emission line pro- 
files (especially in the red-sh ifted electron-scattering wings: 
lHamann fc Koesterke l ll998). (ii) bring theoretical spectral 
energy dis tributions closer to the observed infrared and ra- 
dio fluxes jNueis et alJll998t) . 

There are two broad classes of WR stars: WN, where 
nitrogen lines dominate, and WC (WO), with dominance of 
carbon (oxygen). These two classes may then be further di- 
vided according to spectral ionisation related to wind tem- 
perature. Many of the cooler WC subtypes are known to 
be prolific dust makers. There are two situations in which 
a WC star can produce dust: (1) in a presumably sin- 
gle WC star of late spectral subtype (mainly WC9 and 
some WC8) and (2) in a WC+O binary with a relatively 
long-period, often eccentric, orbit and no restrictions on 
the WC subtype . Currently seven WC stars are known 
to be episodic dust makers with periods or suspected pe- 
riods of several years. We refer here to WR19 (WC4), 
WR48a (WC8), WR70 and WR98a (WC9), and thr ee WC7 
stars WR125, WR137 and WR140 JWilliams Ill995l) . All of 
these are confirmed or suspected binaries with elliptical or- 
bits and massive companions. They serve as unique labo- 
ratories, allowing one to study formation of dust in most 
hosti l e environments |Harr ies et al.||2004fr |Marchenko et all 
1 2001 120031: iMonnier et allll999l l2002t iTuthill et al.lll99a 
Iwilliams 11200 ll) . 

As one of the dust-produci ng prototypes, WR137 
(WC7vd+Q9. \vkn der Hucht l200ll) was frequently observed 
in the last two decades. The initial d iscovery that WR137 
forms dust during large IR outbursts dWilliams et allll985l) 
prompted various investigators to obtain observations in 
1997-1998 during the most recent outburst. One of these 
studies based on near-lR HST/NI CMOS2 images (aroun d 
the predicted periastron passage, iMarchenko et al.lll99S() . 
revealed dust clumps leaving the system. Different IR and 
UBV photometry campaigns yielded a ~ 13 year modulation 
JPanov et alJl200d IwTlliams 11200 lD which was supposed to 
be related to the presence of a binary. WR 137 also shows 
a strong polarisation "line effect", po ssibly produced in a 
spatially anisotropic (flattened) wind jHarries et al.l l2000). 
Two of the best st udied "line effect" WR systems. WR 6 
jHarries et alll998lll999l) and WR134 jHarries et alll998h . 



reveal the p resence of periodically repeatable structures in 
their winds jMorel et al.ll997lll999l and references therein). 
We organised an intense multi-site spectroscopic monitoring 
campaign of WR137 in 1999-2000 in order to improve the 
orbital parameters and study the structure of the WR wind 
and thus obtain a better understanding of the dust forma- 
tion process. Here we report on the results of this campaign 
and its integration into previous observations. 



2 OBSERVATIONS 

Spectroscopic CCD data were obtained from 1986 to 2000 at 
six different observatories (see Table0 with 1-2 m class tele- 
scopes: the Observatoire de Haute Provence (France), here- 
after OHP; the Observatoire du Mont Megantic (Quebec), 
OMM; the Tartu Observatory (Estonia), TO; the Domin- 
ion Astrophysical Observatory (Victoria, Canada), DAO; 
the Ritter Observatory (Toledo, Ohio), RO; and the Kitt 
Peak National Observatory (Arizona), KPNO. 



3 REDUCTION AND ANALYSIS 

All spectra were extracted locally at each observatory us- 
ing standard routines (mainly de-biasing, flat-fielding and 
conversion to a one-dimensional spectrum). The remaining 
calibration steps (wavelength calibration, continuum recti- 
fication, removal of the telluric spectrum) were performed 
for all the spectra in succession using standard 1RAF 1 rou- 
tines, thus producing a homogeneous data set. During the 
continuum rectification we made sure that all the parame- 
ters, especially the wavelength windows assigned to a line- 
free continuum, were identical for all the spectra, although 
the much higher resolution but limited spectral range of the 
KPNO spectra forced us to choose fewer continuum windows 
(although similar ones for the overlapping range) for their 
rectification. Spectral samples used for the determination of 
the continuum are shown in the top panel of Figure 

To be able to detect and measure subtle temporal vari- 
ations in the WR line profiles, we carefully removed the 
telluric lines from all the spectra. Using normalized spec- 
tra of WR 137 taken at different air masses, we created a 
template of the telluric lines, then divided the appropriately 
scaled and binned (to attain the same resolution as the tar- 
get spectrum) template into each observed spectrum. Addi- 
tional iterative adjustment of the intensities and positions of 
the telluric lines allowed us to achieve an optimal removal. 

3.1 Assessing variations in the spectra 

3.1.1 Temporal Variance Spectrum 

Observing line profile variability in WR spectra, one may 
encounter several types of intrinsic variations: a stochas- 
tic componen t related to the small-scale inhomogeneitie s 
in the wind jMoffat et all Il988l: ILepine fc Moffat I ll999T) . 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the Na- 
tional Science Foundation. 
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Table 1. Summary of Observing Runs. 



Site 


telescope 


UT dates 


No. Spectra 


Sp. Coverage 


Sp. Dispersion 


S/N 




(m) 






(A) 


(A/pix) 


(cont.) 


DAO 


1.8 


08/1986-07/1999 


41 


5300-6100 


0.75 


250 


KPNO 


0.9 


11-19/09/1999 


58 


5600-6000 


0.10 


100 


OHP 


1.5 


01/1997-07/1999 


50 


5200-6100 


0.35 


60 


OMM 


1.6 


07/1991-09/1999 


55 


5100-6600 


0.65 


80 


RO 


1.0 


06/1999-09/1999 


67 


5100-6100 


1.50 


90 


TO 


1.5 


06/1999-03/2000 


27 


5600-6000 


0.80 


60 


Total 






236 









as well as periodic variations on various timescales from 
hours to weeks, presumably caus ed either by pulsation, ro - 
tation or an orbiting companion feratschi fc Blecha Ill99(l) . 
To assess whether the variability patterns are intrinsic to 
the star or due to spurious features introduced during the 
reduction process (mainly associated with the continuum 
fitting and wavelength calibration), we use the "Temporal 
Variance Spectrum" (TVS) approach iFullerton et alJll996t 
iPrinia et all Il99fj iMarchenko et alJ Il998i) . At each wave- 
length (index j) we calculate: 



TVS, 



N 



i N - 

— l en 



s 



(Si. 



Si 



-2 2 



(1) 



where index i refers to a given spectrum (total number N), 
(^■) 2 is taken to be 1 (oo is the noise determined for the en- 
tire time series and Oi the noise determined from continuum 
pixels in each spectrum), Sij is a spectrum taken at time 
ti at the j th pixel or wavelength and Sj is the spectrum 
averaged over time (index i). The last term of the equation 
helps to eliminate the spurious details arising from small po- 
sitional errors introduced during the wavelength calibration: 
a is proportional to the typical error of the radial velocities, 
and pj is the gradient of Sj with respect to the wavelength. 

As defined, the TVS does not provide any detailed in- 
formation about temporal characteristics of the line profile 
variability. This can be remedied by calculating TVSs for 
isolated sub-sets. We have therefore calculated a TVS for 
each set of observations. We grouped the data sets in dif- 
ferent ways to access different variability periods. Data sets 
were first grouped on a day-by-day basis, showing significant 
variability of the strong emission lines CIII A5696 and CIV 
A5802/12. Then the data were grouped in gradually larger 
intervals. The most representative sub-set for OHP, OMM, 
DAO, TO data (1999, July) is shown in Figure □ along with 
the TVS for the KPNO data (1999, September). The TVS 
for the RO data is similar in shape but the variations are 
less significant, as expected from the higher noise level of 
the data. It is not shown here for reasons of clarity. 

Statistical behavior of the TVS follows a o"oXiv-i dis- 
tribution, hence thresholds (straight lines in our case) of a 
given statistical significance can be drawn for each TVS: 
where N is the number of spectra and (To 
refers to the noise level in the adjacent continuum. In Fig- 
ure one can see that the variations in the CIII A5696 and 
CIV A5802/12 lines are significant at the 99% level (i.e., ex- 
ceeding the 1% threshold) for the combined OHP, OMM, 
DAO, TO data. 




5400 



6000 



5600 5800 

Wavelength (A) 

Figure 1. The upper panel shows the reference mean spectrum 
of WR137 (combined OHP, OMM, DAO, TO and RO data sets) 
with the superposed, vertically shifted (for clarity) mean spec- 
trum from the KPNO set. The continuum windows used for both 
sets are displayed in the form of horizontal bars. The upper ones 
refer to KPNO data while the lower ones refer to OHP, OMM, 
DAO, TO and RO data. The middle panel shows the TVS for 
the KPNO data for September 11-19 1999, and the lower panel 
shows the TVS for the combined July 1999 sets from DAO, OHP, 
OMM and TO. 



It turns out that the observed profile variations (Fig. 
and other TVSs not shown here) are practically always sig- 
nificant at the 99% level. The spike-like features around the 
sharp interstellar Na D1/D2 lines at A 5890/96 are caused 
by combinations of spectra with different spectral resolu- 
tions (bottom panel) or difficulties in fitting the continuum 
(middle panel). Indeed, the RO TVS (not shown here) does 
not display these features as long as it is a combination of 
data from the same instrument and the continuum is well 
defined. 
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Figure 2. Sequence of residuals of the strongest lines (CM A5696 
and CIV A5802/12) after the average spectrum of the 1999 multi- 
site run has been subtracted from each individual spectrum. The 
residuals are displayed as a greyscale plot and stacked sequentially 
from bottom to top. The white horizontal bands separate spectra 
obtained on different nights. The grey horizontal bands represent 
missing spectra. This Figure represents nights from July 21 to 
31, 1999, where we have the highest density of data. The average 
spectrum is shown in Fig[3] The left panel shows the residuals 
after subtraction, while the right panel shows the same residuals 
after the large-scale variable components have been filtered out 
with the power spectrum shown in Figure [3] Note the variability 
occurring in the form of narrow moving features in the CM A5696 
emission line, and in the form of global line-strength variations in 
CIV A5802/12 on the left panel. While the former is intrinsic to 
the star, the latter is an artefact of the spectral rectification. In 
the right panel, much of the (spurious) global variability in the 
emission line strengths has been removed, leaving only the signa- 
ture of the intrinsic narrow moving subpeaks. Variability patterns 
from the narrow subpeaks are quite obvious in CM A5696 and 
CIV A5802/12. 



3.1.2 Wavelet analysis 

In order to investigate the variability in the different 
spectral lines on a relatively small scal e, it is possibl e to us e 
the wavelet techniques described in iLepine et al] fcOOd) . 
To perform this analysis of rapid spectral variations, a 
carefully computed mean spectrum must be subtracted 
from the individual spectra. 

The mean spectrum was constructed in steps. First, 
the sub-sets of the data (usually the rectified, telluric-free 
spectra from the same observatory) were grouped into 
preliminary means. Then, these groups were cleaned from 
any deviating spectra (comparing the preliminary mean 
to the individual spectra), forming the final mean spectra. 
The final mean spectra at similar spectral resolution 
(OHP, OMM, DAO, TO) were used to form, via simple 
averaging, the template mean used in this section. Due 



Figure 3. Mean wavelet power spectrum of the residuals (top) 
and of the wavelet filtered residuals (centre) . The residuals yield 
a strong response at small scales near the CM A5696 line, and 
at larger scales near CIV A5802/12, consistent with the variabil- 
ity patterns observed in Figure [5] Filtering out the large-scale 
(instrumental) components leaves only the signature of the (in- 
trinsic) narrow variable subpeaks. The wavelet power spectrum 
of the filtered residuals yields a positive identification of narrow 
emission features in both CM A5696 and CIV A5802/12, and sug- 
gests that narrow variable features are also present in Hell A5411, 
CIV A5471 and OV A5592 (grey patches). 

to the different spectral resolutions, the KPNO and RO 
data were first treated separately. We show that all the 
most prominent emission lines in the 5400-6000A range 
are variable, predominantly in the form of narrow emission 
subpeaks moving across the broad line profiles. 

In Figure (left) we plot deviations from the average 
profile (shown at the bottom of Figure |3J for spectra ob- 
tained during the 1999 multi-site campaign (OHP, OMM, 
DAO, TO, RO) for the CIII A5696 and CIV A5802/12 lines 
and only for the nights where the temporal coverage is the 
best. Two types of variability patterns are apparent in the 
residuals: (1) patterns marked by narrow moving features 
as is most apparent in the CIII A5696 emission line, and 
(2) broad patterns that suggest that the whole line is glob- 
ally increasing or decreasing in strength, as is most obvious 
in the CIV AA5802/12 emission doublet. The first patterns 
correspond to the moving emission-line features s een in the 
spectra of many WR stars iLepine fc Moffat 119991) . The sec- 
ond pattern, on the other hand, is largely an artefact of the 
data reduction procedure. All spectra have been rectified to 
a pseudo-continuum level. The pseudo-continuum is poorly 
defined in the regions with numerous blended, broad emis- 
sion features. The chosen pseudo-continuum sampling pro- 



Spectroscopic Study of HD192641 5 




5650 5700 5750 5800 5850 5900 




5650 5700 5750 5800 5850 5900 5950 
wovelength (A) 



Figure 4. The KPNO run: temporal sequence of residuals af- 
ter the large-scale variable components have been filtered out - 
compare to Figure 121 



duces relatively small, but, nevertheless, easily detectable 
systematic differences between the groups of spectra (usu- 
ally coming from different observatories). 

Fortunately, the variability of the intrinsic narrow mov- 
ing features and that of the global (spurious) line strength 
variations have very different characteristic scales. It is thus 
possible to separate them using wavelet analysis. In Figure|3J 
we display the mean wavelet power spectrum of the residu- 
als. _Th^jwrowJ2isCTete_Wind Emission Elements (DWEEs; 
see iLepine fc Moffat Ill999l) introduce a significant response 
on scales of ~ 1 A. Note the very dark patch near 5700A 
that is associated with the high-amplitude DWEEs in CIII 
A5696. Global line variations, on the other hand, yield re- 
sponses over much larger scales. The very large, dark patch 
near 5800A that extends to scales of several tens of A is the 
signature of the large pseudo-continuum rectification arte- 
fact in CIV AA5802/12. Note also the narrow dark patch 
near 5890 A extending over a range of very small scales (< 1 
A): this corresponds to another artefact associated with the 
interstellar absorption doublet. Because the spectra from 
different observatories have slightly different spectral res- 
olutions, subtraction of the mean profile from the individual 
spectra introduces residual features that are very apparent 
in Figure [5] 

We use wavelet filtering to remove the rectification arte- 
fact from the time series of residuals. Keeping only wavelet 
coefficients on the smallest scales (< 4A), we use the wavelet 
reconstruction theorem to rebuild each individual residual 
spectrum. The mean wavelet power spectrum of the filtered 
residuals is displayed in Figure^] A time series of the filtered 
residuals is displayed in the right panel of Figure H and can 



be compared to the original time series in the left panel of 
the same figure. One sees that all large scale variable features 
have been filtered out. The narrow moving features in the 
CIV AA5802/12 emission doublet, which initially were com- 
pletely drowned out in the (spurious) global line-strength 
variations are now clearly visible. These intrinsic features 
are weaker than those in CIII A5696. This behavior mimics 
the variability pat terns seen in emiss ion lines of the single 
WC8 star WR135 jLepine et al.ll200rjT) . 

A thorough examination of the filtered residuals and 
corresponding wavelet power spectrum also reveals what ap- 
pears to be a very weak signature from intrinsic narrow emis- 
sion features in the three prominent emission lines bluewards 
of CIII A5696. In Figure^! darker patches can be seen match- 
ing the locations of Hell A5411, CIV A5471 and OV A5592, 
indicating a significant component of variable features at the 
scale expected for the narrow emission-line subpeaks. 

The same wavelet-filtering procedure is repeated for 
the spectroscopic time-series obtained during the Kitt-Peak 
run. Again, random variations on large scales affect the CIV 
AA5802/12 doublet, and these are an artefact of the spectral 
rectification. Fortunately, these artefacts occur on a scale 
(> 10A) that is larger than that of the intrinsic DWEEs, 
and they are easily filtered out. The resolution of the Kitt- 
Peak spectra (0.1 A per pixel) is higher than that of the 
other sub-sets of data, hence one may obtain a clearer pic- 
ture of the scale properties o f the DWEEs. The behav ior is 
largely the same as noticed bv lLepirie fc Moffat Hl999l) . pro- 
ducing, in the wavelet power diagram, the shape of a sym- 
metric triangle standing on its tip which marks the centre of 
the emission line, where the DWEEs are narrower than near 
the edges of the line. The time-series of the wavelet-filtered 
KPNO spectra is shown in Figure 01 The general motion 
of the DWEEs, from the line centre toward the line edges, 
is quite apparent in CIII A5696 and CIV AA5802/12, thus 
confirming the presence of a strong variability as seen in the 
TVS analysis in section 13.1.1 1 From the KPNO time-series, 
as for the previous one, the DWEEs in CIII A5696, if com- 
pared to CIV AA5802/12, do not give the impression that 
they follow similar patterns. However, the doublet nature of 
the CIV transition may confuse things considerably. Hence, 
we must explore the possibility that the variability in the dif- 
ferent emission lines follows similar patterns. Such behavior 
would be consistent with a model of the wind where the dif- 
ferent emission lines are formed at different depths in the 
wind on average, depending on their excitation potential, 
but with a significa nt overlap between the emission regions 
jHerald et alj200Cl) . Therefore, line profile variations caused 
by outmoving clumps wo uld be expected to show up in a ll 
the lines at similar times iLepine Ill996l : ILepine et alj f2000). 
This will be discussed in detail in section 13.2.21 

3.1.3 Periodograms 

Inspired by the positive results of the TVS analysis, we 
also searched the spec t ra for periods using the periodogram 
approach of IScarglc (1982) combined with the cleaning 
algorit hm of iRoberts et all (I1987T) . Following IScargle I 
(1982), we use iVo = N/2 frequencies uj, noting that the 
number of observations, N, could be different for the TVS 
and periodogram samples. This defines the corresponding 
false alarm probability thresholds: T ~ ln(No/po), where 
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Figure 5. Kitt Peak Observatory periodograms reduced to one 
dimension. The average periodograms for the P Cygni absorption 
part of the three main lines (CHI at 5670-5692, CIV at 5755- 
5776 and Hel at 5838-5848) are shown in the left panel while 
the average periodograms of the two main lines (CIII at 5670- 
5720 and CIV at 5785-5835) along with the continuum (AA 5645- 
5665,5735-5755,5900-5930) are shown in the right panel. The false 
alarm probability threshold levels are also indicated. 



po is the false-alarm probability, usually set at 0.01-0.05. 
The cleaning procedure effectively suppresses the false 
periodicities arising from the non-regular character of 
the observations, the most prominent ones coming from 
the daily sampling bias of the data sets (see Fig. |SJ. 
Although we found no periodic behavior in the spectra 
from OHP, DAO, TO, RO or OMM, this method yielded 
a small-amplitude, 0.83 ± 0.04-day period in the KPNO 
spectra, searching in the range P=0.2-30d. Probably, both 
the high homogeneity of this particular set, as well as its 
high spectral resolution helped to detect the low-amplitude 
signal. 

We use wavelet filtering to isolate two different compo- 
nents of the variability: small-scale 4 A) and large-scale 
(> 4 A) structures (methods described in section 111 1.211 in 
order to distinguish between real variability and artefacts 
of the spectral rectification. The presence of this 0.83d 
period in the large-scale structures (i.e. associated with 
overall variations of the continuum caused by errors in its 
fit) would invalidate its existence while its appearance in 
the small-scale structures would suggest that it is real. 

We found no significant signals in the wavelet filtered 
(either small- or large-scale) periodograms calculated from 
averaged fluxes of the emission parts of the three major lines: 
CIII (5670-5720 domain), CIV (5785-5835) and Hel (5865- 
5885). On the other hand, the periodogram of the P Cygni 
absorption parts shows a clear peak at 0.83d for CIII while 
the peaks are diminished for the CIV and Hel lines. More- 
over, the 0.83d signal is absent from the large scale structures 
in the CIV, CIII and Hel lines of the same periodogram. 
Combining the results from Fig|^l and the wavelet filtered 
periodograms, we conclude that: (1) there is no periodic sig- 
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Figure 6. Absorption troughs of CIV AA5802/12 (left) and Hel 
A5876 (right) phased (P=0.83d, T =HJD 2480198) with the 0.83d 
period found exclusively in the KPNO data. Different lines depict 
average profiles calculated for the given phase bins. 
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Figure 7. Averaged flux in the 5750-5780 A range of the CIV 
AA5802/12 line versus phase (P=0.83d; T =2450198d) for the 
KPNO data. 



nal in the emission parts as suggested by the wavelet filtered 
periodograms of the lines; (2) there is a periodic signal at 
0.83d, but it can be seen only in the small-scale structures 
in CIII A5696 (it is amplified by the filtering) and in the 
absorption part of Hel A5876 (Fig. [S] left panel). Note that, 
though the unfiltered signal in CIII A5696 falls much below 
statistically significant levels, wavelet filtering allows one to 
increase the period's visibility rather dramatically. This dif- 
ferent behaviour of the CIII line is to be expected as it has 
no discernible P Cygni trough, and all the variability seems 
to be concentrated there. It is no surprise either to find no 
0.83d signal in the CIII line in the unfiltered data: this signal 
shows only in the filtered small-scale set, being restricted to 
the blue wing of the CIII line profile, roughly located at the 
same place as the CIV and Hell counterparts. 

We show the phase-binned profiles of CIV and Hel in 
FigHJand an average flux coming from the P Cygni absorp- 
tion trough of CIV in Fig. |7| Both the amplitude (0.5%- 
1.5% of the continuum) and period of this variation are 
small compared to similar cases as seen for two other sin- 
gle WR stars with similar large-scale periodic spectral varia- 
tions (WR134 and WR6), where the amplitude is 5%-10% of 
the co ntinuum and the periods are 2.31d and 3.76d, respec- 
tively (iMorel et alJll997L 11999^ . Note that, by themselves, 
the significance levels plotted in Figure|5]are uncomfortably 
low. This may be related to the relatively high level of the 
noise introduced by the adopted procedure of rectification: 
note that, while calculating the respective values, we used 
the unfiltered KPNO spectra. However, Figures [(J and [7| lend 
strong support to the 0.83-day periodicity. Note that these 
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Figure 8. Clump evolution during two nights: left panel - 
21/07/1999; right panel - 22/07/1999. The mean profile for all 
runs for the CIII A5696 line is shown in the upper panels. The 
bottom panels show montages of the deviations from the mean 
profile. In this plot we trace only the longest-living structures 
(straight dashed lines). However, practically all moving features 
were chosen for further analysis (Fie. 19*1, 



periodic variations show coherent patterns in two regions 
separated by ~ 50A (~ 250 resolution elements). We fur- 
ther discuss the 0.83± 0.04d period in section FOl 



3.2 Spectral-line variability 

3.2.1 The (3 wind-velocity law 

We then examined the rapid spectral variations in order to 
try and constrain the f3 wind-velocity law. To achieve this 
goal, the mean spectrum mentioned in section 13.1.21 was 
subtracted from the individual (unfiltered) spectra. Figure 
|H] displays examples of the spectral variability of the CIII 
A5696 emission line for the nights of July 21 and July 22, 
1999. With the removal of the mean profile (top panels), 
it is much easier to follow relatively weak emission details 
travelling on top of the broadest lines (for similar results 
in other stars see iGrosdidier et alJl200lt iLepine fc Moffat I 
1999). As can be seen in Figure |5] individual subpeaks can 
be followed typically for several hours. It is customary to 
relate the small emission peaks to small-scale density en- 
hancements ("clumps" or "blobs") moving with the WR 
wind. Blobs are assumed to follow a f3 wind- velocity law 
like the rest of the wind (in fact they may constitute the 
wind as a whole): 

v(r) = Voo (l-—f. (2) 



Their acceleration is 

The latter is equivalent for (3 > 2 to iLepine fc Moffat 119991) : 

a(v) = [vln(— )f/PR*- (4) 

Assuming = 1885 kms' 1 and R+ = 4.5 Rq for WR137 
(see Tabled one can obtain an estimate of (3. 

Following individual subpeaks (Fig. [HJl, one can draw 
straight lines which have a slope proportional to the ra- 
dial acceleration (i.e. the actual acceleration projected by 
an angle 8 onto the line of sight) of the moving features. 
Thus, measured accelerations (or) along with the derived 
radial velocities (vr) are plotted in an acceleration-velocity 
diagram (Fig. |UJ, which is traditionally used to constrain 
parameters of the j3 wind-velocity law. The RVs of the sub- 
peaks have been measured by fitting gaussian profiles to the 
moving peaks, assuming a(RV) ~ 50 fcms -1 for the blobs 
in each spectrum. The least-square fitting of a straight line 
to the positions of the peaks provides estimates of the stan- 
dard deviations for the acceleration. Due to the exceptional 
temporal coverage of the two nights represented in Figure 
|H| it was possible to follow subpeaks even in the regions cor- 
responding to the steep flanks of the CIII A5696 line. These 
features move with greater radial acceleration than the blobs 
originating near the centre of the line, due to projection ef- 
fects, i.e. one observes an = a(v)cos8 and Vr = v(r)cos9 
at #=const. for each subpeak, such that those at line centre 
have 8 ~ 90°, and those near line edge have 8 — > 0° on the 
blue side or 8 — * 180° on the red side. 

Figure [5] shows that most of the components satisfy the 
expected clrVr ^ relation, thus confirming the initial as- 
sumption that the moving peaks are related to radially accel- 
erated structures leaving the star. Note that measurements 
with excessively large errors in acceleration, a aR ^ 0.02 
kms~ 2 , have not been drawn in Figure|U]to avoid clutter. 
The 142 spectra (10 nights of observation at OHP, OMM, 
DAO, TO and RO during the period of 1999 July 21-31, 
except for July 29 when only one spectrum was available, 
plus 8 nights of obervation at KPNO during the period of 
September 1999) result in /3R* = 25 ± 5 R Q , i.e. j3 ~ 5, 
which is reaso nably close to the r ange of values determined 
for WR137 bv lLepine fc Moffat I <1999h {fiR* ~ 35-90 R Q ). 

Note that the lower-acceleration structures outnumber 
the higher-acceleration ones. One may find an explanation 
for this in Figure [TOl which shows the actual density of blobs 
compared to their theoretical density (N(Q) ~ sin(8), as- 
suming constant spherically-symmetric density) versus the 
angle with respect to the line of sight (8). The theoretical 
density maximum has been adjusted by a multiplicative fac- 
tor to the maximum of the observed distribution. As can be 
seen, both observed and theoretical densities fall off as one 
goes away from #=90°. However, there seems to be a real 
deficit in the observed density of blobs coming toward us (8 
~ 0°) or going away from us (8 ~ 180°) compared to the ex- 
pected density. This can be explained by the fact that events 
responsible for the high acceleration structures normally oc- 
cur far from the centre of the spectral line, thus being dif- 
ficult to observe and track: (1) the corresponding features 
move rapidly and require good temporal coverage and (2) 
they are hard to recognize and trace because of their po- 
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Figure 9. Radial acceleration versus radial velocity for detected blob trajectories from the DAO, OMM, OHP, RO, TO spectra of CIII 
5696 for the period of July, 1999. The horizontal bars span the range of observed RVs while the vertical bars represent ± the standard 
deviations for the fitted accelerations. A few points, due to measurement errors, fall into "forbidden" regions, i.e. on the back-side of the 
wind coming toward us and vice versa, providing wrong-sign accelerations. Curves based on eq. (4) represent angles to the line of sight 
from 8=0 to 180 degrees in steps of 30 degrees, left to right, bottom to top. The dash-dotted lines represent the inner and outer limits of 
the le r in velocity (V e = 1180 ± 10 km s~ 1 and AV e = 270 ± 10 km s~ 1 for the CIII A5696 line) calculated according to lLepine fc Moffat I 
( 1999) for (3R* ~ 25Rq. The two zones delimiting the ler in velocity space are indicated by horizontal bars. 



Table 2. Stellar parameters of HD192641. 



distance 


R* 


P 




Mwr 


M 


Mwr/Mo 


P L X 


(kpc) 






(km s —1 ) 


(Ma) 


(Ms) 




years (erg s~ ) 


1.82" 


4.5 ± 2.5 a 


5.6 ± 1.2 6 


1885 ° 


4.4 ± 1.5 6 


20.0 ± 2.0 C 


0.22 ± 0.07 6 


13.05 ± 0.18 6 0.25 ± 0.14 10 32d 



(aJ_^^^^^^^^g3J^M9^^ggteri^^^^fc^^^^3 > v °° c l° se to 1900 kms 1 from fvan der Hucht I 1200 ll) : (b) this paper: (c) 

IVacca et al] <1996t)lGies I i2003h : fdl iPollock etldTll995ft . 



sitions and broader/flatter profiles superposed on the steep 
slopes of the line. In Figure [5] the limits of the CIII A5696 
line emission region (ler ) have been calculated according to 
iLepine fc Moffatl 1^999). Note that some points are located 
well outside the central cluster of points in Figure |0J i.e. 
outside the ler in the "allowed" regions. These "outliers" 
might result from the superior temporal coverage provided 
by the more complete data from July 21 and 22 and some 
nights from the KPNO dataset. This may have allowed one 
to follow weak blobs on the steep flanks of the lines, with 



stochastic errors putting the odd point outside the allowed 
area. 



3.2.2 Cross-correlation Analysis 

We now verify if the variable subpeak patterns in the differ- 
ent emission lines are similar. A visual examination of the 
spectroscopic time series (Figs. [5] does suggest that the 
variability pattern in different lines may indeed be similar, 
i.e. that specific narrow subpeaks show up in each line, at the 
same position of the line profile, and at any given time. To 
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Figure 10. Actual density of blobs (histogram) compared to the- 
oretical density (continuous line) versus angle of projection. See 
text for details. The relative deficit of observed blobs vs. theoret- 
ical number of fast structures towards the wings is likely due to 
the difficulty of measuring blobs in the wings of the CIII A5696 
line. 
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Figure 11. Results from the cross-correlation analysis of the vari- 
ability patterns observed in the CIII A5696 emission line. The left 
panel shows the 1999 multi-site campaign (OHP, OMM, DAO, 
TO, RO) results while the right panel shows the results obtained 
for the KPNO run. Local maxima indicate regions where the vari- 
ability patterns are correlated to those detected in CIII A5696 
(note the value of 1.0 reached at 5696 A indicating maximum 
correlation of the pattern with itself). Vertical lines mark posi- 
tions of the most prominent emission lines relative to the part of 
CIII A5696 used for the correlation. Upper panels show the cor- 
relation functions calculated using the complete pattern in the 
CIII A5696 line (5676< A <5716). The maxima are higher when 
one uses the patterns from the red side (5696 < A < 5716) of the 
CIII line (bottom correlation panels). The correlation essentially 
vanishes when one uses the blue-side patterns, 5676< A <5696 
(middle correlation panels). 



quantify the similarity in the subpeak patter ns, we apply 
the cr oss-correlation technique developed by ILepine et alJ 
(2000). The well-defined variable subpeak pattern in the 
CIII A5696 line is used as a reference (the correlation re- 
sponse). Properly normalized, perfect match of the pattern 
at all times would yield a correlation value of 1.0, while the 
absence of any similarity should yield a correlation value 
close to 0.0, and —1.0 for an anti-correlation. Results are in- 
dependent of the intensity of the pattern; hence if a pattern 
is identical to the reference pattern, it will yield a high corre- 
lation value, even if its amplitude is only a fraction of that of 
the reference pattern. One caveat is that instrumental noise 
can significantly degrade a subpeak pattern, especially if its 
intensity is low. Since the variable subpeak patterns in Hell 
A5411, CIV A5471 and OV A5592 are near the noise limit, 
one should not expect them to be highly correlated with 
CIII A5696. However, a small, but clear correlation should 
be regarded as significant. 



In order to minimize the degradation in the correlation 
because of instrumental noise, we apply the cross-correlation 
technique to the wavelet-filtered spectroscopic time series 
(Figs. Results of the cross-correlation analysis are dis- 
played in Figure ITT1 The cross-correlation is performed sep- 
arately on the KPNO time series, compared to the multi-site 
(OHP OMM, DAO, TO, RO) time series. As in lLepine etafl 
(2000), we use three different parts of the CIII A5696 line 
as the reference. First, we use the complete pattern in the 
CIII A5696 line (5676< A <5716) as our correlation response 
(Fig. 1111 top). As expected, we find a correlation value of 
1.0 at the position of the CIII A5696 line itself. We also 
find secondary maxima at the exact positions of all the most 
prominent emission lines. Although some of these maxima 
are quite weak (with a correlation value of « 0.1) they ap- 
pear to be statistically significant because they occur at the 
exact wavelength positions expected for these atomic lines. 
The low correlation coefficient is a direct consequence of the 
weakness (proximity to the noise level) of the variability pat- 
terns in all the lines, except in CIII A5696. It is also possible 
that the variability patterns in the different lines are not 
exactly the same, as, generally, the lines do not form in the 
same parts of the wind. 



We repeat the cross-correlation analysis, this time us- 
ing only half of the variability pattern, either on the blue 
(5676< A <5696) or red (5696< A <5716) side of the CIII 
A5696 emission line (Fig. 1111 second and third panels from 
the top). Results show a much stronger correlation when one 
uses the red side of the emission line. This perfectly mim- 
ics the behav ior of the subpeak s in the presumably single 
star WR135 (ILepine et alJl2000l) . Most WR emission lines 
are flanked on their blue side by a weak but variable P- 
Cygni absorption trough, with a notable exception of the 
CIII A5696 emission line. Generally, the P-Cygni absorption 
trough is the most variable part of the profile. This may in- 
troduce an extra component of variability on the blue edge 
of the line. However, because CIII A5696 is unaffected by 
P-Cygni variations, the intrinsic variability of its blue edge 
is significantly different from that in the other lines. This 
explains the lower correlation between the blue edge of CIII 
A5696 and the blue edge of the other lines. 
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Figure 12. Radial velocities of the WR component pha sed with 
the 4766d o r bit are shown in the u pper panel: data f r om | A nnuk 
Jl99ll.ll995h . lMassev eta!] <198ll) and iMoffat et all Jl986ft refer 
to photographic plates; the rest are measured from the current 
CCD spectra from DAO, OMM, OHP, RO, TO and KPNO. A 
typical standard deviation for the measurements is ~ 15 kms~ 1 
for the former and 10 kms~ 1 for the latter. The bottom panel 
shows radial velocities for the O component. Full lines show the 
derived orbit. 



3.3 Radial velocity variations: Orbit of the two 
components 

Radial velocities (RV) were measured by fitting gaussian 
profiles to the Hel A5876 absorption line corresponding to 
the O star and cross-correlating the template (an average 
constructed from the data sets coming from different obser- 
vatories) with the individual spectra in the region covering 
the strong CIV AA5802/12 emission blend for the WR com- 
ponent. Typical uncertainty of the measurements run at ~ 5 
- 10 km s" 1 for an individual spectrum of the WR star and ~ 
10 - 20 kms~ x for the O component. The results were cor- 
rected for systematic differences between observation sets. 
This adjustment was done taking the RO and OMM ob- 
servations together, as in the overlapping epochs they show 
negligible relative shifts for the WR star measurements: the 
average velocity measured on the overlapping sub-sets (sum- 
mer 1999) was < V RO > = 3.0 ± 6.7 kms' 1 , < Vomai > = 
3.6 ± 5.6 kms" 1 . Then, we shifted the Kitt peak RV data 
into the RO and OMM system and applied the same pro- 
cedure to the less numerous DAO and OHP sets (shifted 
relative to the RO+OMM system by 7 - 9 fcras" 1 ). The 
Tartu set was not modified, as many of its data were taken 
later on, without any substantial overlap (see Tabled pub- 
lished separately in the electronic version). 

With the assumption that our data set alone may 
not be adequate in revealing the presumably long orbital 
period of WR137, we have incorporated earlier published 
RV measurements (iMassev et al.lll98ll : lAnnuk lfl99ll Il995t 
Moffat et al]| l986|). We applied t he orbit-fitting algorithm 
from iBertiaufe Grobben"l 1^969), assigning individual 7- 
velocities for the following sets: all the new CCD data for the 
WR component, as discussed above (made ho mogeneous by 
additional adjustments), then the WR sets of lMassev et alJ 



Table 4. Orbital elements of HD192641 


Parameter 


WR O 


P(d) 


4766 ± 66 


T (JD) 


2450198 ± 186 


*> (°) 


326 ± 15 


e 


0.178 ± 0.042 


7 (km s -1 ) 


-14.6 ±0.9 -5.4 ±0.7 


K (kms- 1 ) 


27.9 ±1.2 6.1 ±1.3 


°CCd{o - c) (kms- 1 ) 


6 15 


^plates (0 - c) (kms' 1 ) 


11 26 


asini (AU) 


12 ±0.5 2.7 ±0.6 


Msin 3 i (Mq) 


3.4 ±1.0 15.3 ±2.1 



These parameters were derived with the combined radial veloc- 
ities of CIV A5806.0 and absorption at H e l A58 75.6 using the 
program developed bv lBertiau fc Grobben I Il969fl . 

l|l98ll) and lMoffat et alJ jl98ff l (photographic data, treated 
as 2 individual sub-sets) , then 2 sets for the O star, one from 
our ca mpaign, and the second described in IMassev et alJ 
(1981). This allowed us to compensate for the systematic 
differences in the RVs coming from different systems and 
epochs and evaluate the overall accuracy of the data: for 
the CCD data, <t _ c is ~ 6 kms -1 for the orbital fit of the 
WR component while it is ~ 15 kms -1 for the O compo- 
nent, in line with the accuracies attained in measurements of 
individual spectra. The photographic plates yielded <r _ c ~ 
11 km s _1 and ~ 26 kms -1 for the WR and O components, 
respectively. Despite the different accuracy of the photo- 
graphic and CCD data, we have kept the same weights while 
deriving orbital elements: less accurate photographic data 
have, nevertheless, an advantage of being taken at differ- 
ent epochs, while the CCD set mainly concentrates around 
1998-2000, which is not very helpful in the situation when 
the orbital period exceeds 10 yrs. 

The orbital elements obtained from the least-square fit 
of 427 WR RVs and 373 O RVs are listed in Table H Figure 
H2l shows that the computed orbit based on these parameters 
fits quite well all RV data obtained so far for the WR compo- 
nent. On the other hand, the fit for the O component is less 
precise, due to relatively large errors of the measurements. 



4 DISCUSSION AND CONCLUSIONS 
4.1 Spectral line variability 

The presence of propagating stochastic structures indicates 
a clumpy wind as in all WR stars appropriately studied so 
far. Moreover, this study enabled us to constrain the param- 
eters of the assumed /3-law with good accuracy, due to the 
large quantity of data. We found /3i?* ~ 25 ± 5Rq (accord- 
ing to Figure [5J which agrees marginally well with previ- 
ous results, ~ 35-90 Rq llLepine fc Moffati ri999). found us- 
ing a wavelet analysis for WR137. However, one should be 
cautio ned against any overs implified interpretation of the 
/3-law. lHamann et al.l <l200ll) have discussed the apparent 
contradiction between the fast velocity law ((3 ~ 1) and 
the slo w drift of the DACs (Disc rete Absorption Compo- 
nents: ICranmer fc Owocki Ill996i) in the winds of O stars. 
The explanation resides either in a difference between the 
velocity of patterns (DACs) and structures (modulations) 
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Table 3. Radial velocities of the WR and O components of WR137. 



HJD- 


RV 


HJD- 


RV 


HJD- 


RV 


2440000 


km s~ 1 


2440000 


km s -1 


2440000 


km s~ 1 


WR 


Comp. 






<r~ 5-10 km s~ 1 




KPNO 


n=58 


OHP 


n=50 


RO 


n=67 


11432.63091 


0.2 


10836.23647 


21 


11337.86891 


-9.2 


11432.66730 


1.1 


10837.24652 


21 


11350.80538 


-5.1 


11432.74425 


3.5 


10461.23148 


15.54 


11355.78767 


-13.4 


11432.76590 


4.7 


10454.24317 


11.9 


11364.71949 


0.6 


11432.78742 


5.6 


10456.22850 


11.15 


11364.84482 


-11.1 



Table 151 is published in its entirety in the electronic edition of this article. A portion 
is shown here for guidance regarding its form and content, n corresponds to the total 
number of spectra for one observatory. 



in the wind or in different velocity laws co-existing at the 
same time in different spatial zones of the wind. Indeed, 
one may consider alte rnatives to the original /3-law. E.g., 
iHillier fc Miller] il999T) introduced two power indices, thus 
dividing the wind into two distinct parts with different ex- 
pansion laws. Although such a law might fit real data better, 
it is really no more physically justified than the simple /3-law. 

In addition, wavelet analysis and cross-correlation yield 
other interesting results. First, although weak, there is a cor- 
relation in variability patterns between all the most promi- 
nent lines in the AA5300 - 6000 A range. Moreover, this 
correlation is stronger for the red-side than the blue-side 
variability patterns of CIII A5696: in the region of in- 
terest most emission lines -except the strong CIII A5696 
line- show variable P-Cygni absorption troughs, which ef- 
ficiently decorrelates the observed patterns. This behaviour 
was also detected in spectra of the single WC8 star WR135 
jLepine et al.ll2000D . 



4.2 Short period 

The 0.83d±0.04d period could be related to large-scale, co- 
herent structuring of the WR wind, taking the form of 
CIRs (Corotating Interaction Regions: ICranmer fc Owocki I 
Il996l) . CIRs are thought to be responsible for the pe- 
riodic, P=^766d J _rjatterns re-appearing in the wind of 
WR6 JSt-Louis et al.lll995l . their Fig. 5), with a coherency 
timescale of at least ~ 4 rotation cycles. 

If this 0.83d period was related directly to a phe- 
nomenon occuring on the surface of the star, it could imply 
a rotational velocity that is closer to the break-up velocity 
in WR137 than in WR6 or WR134. Taking parameters from 
Tabled we obtain a rotational velocity of 275 kms~ and a 
break-up velocity of 475 fcms" 1 , which give V ro t/Vbreak-up 
~ 0.58. The results are not inconsistent with a rotational pe- 
riod, although one might also want to look for other expla- 
nations. This period could be only a fraction of the real ro- 
tatio n period, as is often supposed for CIRs iSt-Louis et alJ 
Il995l) with multiple arms. Moreover, relatively rapid rota- 
tion of WR137 could be related to a fl attened wind, as im - 
plied by the depolarized emission lines (lHarries et al.l feoOO'). 

The O star companion with its period of 4766 days can- 
not possibly account for such short-term variability, as the 
0.83d period is seen exclusively in the absorption troughs 



of the lines originating in the WR wind. The existence of 
a much closer, low-mass compact companion could be an- 
other explanation. However, WR137 cannot be considered as 
a bri ght X-ray source: Lx (0.2-2.4fceV)=0.25± 0.14 10 32 erg 
s _1 dPollock et al.lll99rj) . On the other hand, the absence of 
a high X-ray output (compared to Cyg X-3 at ~ 10 38 ergs' 1 ) 
does no t necessarily mean that there is no compact com- 
panion fergma fc Yungelson lll998l iMarchenko et af]ll996T) 
in so far as the propeller effect created by the rotating mag- 
netic field may inhibit production of X-rays. A period of 0.83 
days would imply an orbital separation ro — 7 Rq (6.7±0.8 
Rq) between the WR star and a hypothetical companion, 
according to the mass for the WR star (see Table and 
considering a neutron star of 1.4 Mq. Since i?* ~ 4.5 Rq, 
this hypothesis remains plausible until further evidence can 
be found to the contrary. 

Another possible source of the 0.83d variability could be re- 
lated to pulsations of the WR star. Periods of radial pul- 
sations in WR star s are predicted to range from several 
minutes to an hour iBratschi fc Blecha II1996T) and thus can 
be ruled out. The same holds true for strange-mode pul- 
sations iSaio et al.lll998h . for which expected pe riods are 
more likely to be in the r ange of several minutes iSchaller I 
Il992t iGlatzel et aljll993T) . The cla im of P ~ 627s fo r the 
southern Wolf-Rayet star W R40 felecha et all Il992ft was 
contradicted by the re sults of IMarchenko et alJ 1119941) and 
ISchneider et alJ lll994T) . who claimed this period could be 
either spurious (caused by atmospheric transparency varia- 
tions) or a pulsation appearing and vanishing on timescales 
of hours to a day (although no build-up or fading was noticed 
in the intensity of the variation). On the other hand, non- 
radial pulsations with a period ~ 20 hours (0.83d) are a pos- 
sibility. They could be driven by the Edding ton e-mechanism 
jMaeder fc Mevnet l2o"()2llEddington Il9ish . The timescales 
of thes e modes would be of the order of hours to l ess than 
a day jNoels fc Scuflairdl 19861 lDziembowski"lll994ll so that 
they could be held accountable for the observed period in 
WR137. These are related to g+ modes of non radial pulsa- 
tion, but are not believed to last more than a few 10 3 years. 
Further studies of WR137 with good temporal and spec- 
tral resolution (comparable to the current KPNO spectra) 
and high signal to noise ratio (S/N ~ 100-200) could help 
to confirm and determine the origin of this small amplitude 
periodic variability. Note that the non-radial g-mode oscil- 
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lations should have a variable period fealona et al] ll989). 
while the rotationally-modulated variability should provide 
a stable periodic modulation, thus providing means to dis- 
tinguish between the two. 

4.3 Orbital period and masses of the components 

We have computed a reliable orbit for the WR star and 
provided the best possible estimate for the O-component 
orbit. Our derived period of 4766 ± 66 days is in perfect 
agreement with that found independently from t he near-IR 
photo metric data taken between 1985 and 1999 iWilliams I 
l200lft . which yielded 4765 ± 50 days. The small eccentric- 
ity (e = 0.178) seems to be compatible with the ep isodic 
dust formation depicted in iMarchenko et all (^999), who 
assumed e ^ 0.3, based on the rather slow variability of 
the IR flux i Williams 1 1200 ll) and behavior of the ejected 
dust cloud. Assuming that we see the system fairly near to 
edge-on, this would imply that the dust clumps observed in 
1997 September - 1998 May with HST/NICMOS2 are in- 
deed formed around the time when the O-companion was 
relatively close to the WR star. Indeed, the periastron pas- 
sage seems to occur around 1996 (between February and 
October), according to Table [I] Assuming a distance to the 
system of 1.82 k pc (Tabled and a spati al ejection velocity of 
~ 2000 km s _1 iMarchenko et al.ll999ft . the displacement of 
the persistent dust cloud by 0.25" in May 1998 implies that 
the cloud was formed during periastron passage in 1996, 
long before the 1997-1998 images were taken. The dimin- 
ishing orbital separation around periastron should intensify 
the wind-wind collision and, finally, increase the dust out- 
put (e.g. like another long-pe riod dust-producing binary WR 
140: IMarchenko et alJ (|2003)). One more encouraging detail 
is the orientation of the WR orbit. Indeed, the behavior of 
the dust cloud implies that it was ejected from the system 
practically in the plane of the sky , at adust ~ 27° relative to 
the plane iMarchenko e£a0[l999). This fits well the longti- 
tude of periastron passage of the O star: 326° — 180° = 146° 
(Table |3 which implies a orbit = 180° - 146 = 34°, thus 

&orbit ^ Odust^^^^^^^^^^ 

1 ■ | . : ' . : :• - 1 : ' ^ ~1 ||2001) catalogue provides an abso- 
lute magnitude of Mv=-5.7 for the WC7pd+09 system, and 
Mv=-4.52 for the WR-component, im plying Mv=- 5 .25 fo r 
the O-component. Thus, according to IVacca et alJ dl996^> . 
the 09 component should have a spec tral class be tween III 
and V and a mass of Af o =20±2 M Q iGies 112003ft . This, in 
turn, implies Mwr=4A±1.5 Mp, and an incl i nation angle 
of i=67°. On the other hand. lNueis fc Lamersl i2000ft give a 
mean of Mwr=^A ± 1.3 M Q for WC7 stars which, if as- 
sumed to be applicable to the WC7 star in WR137, leads to 
M o =60.9 ± 15.5 M Q and 1=39°. However, the shape of the 
ejected dust cloud strongly favors the larger value of inclina- 
tion. Moreover, the mutual orientation of the orbital plane 
and t he plane containing the flattened wind iHarries et alJ 
12000ft also points to a rather high value of i. 
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Table 3: Radial velocities of the WR and components of WR137. 
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Tabic 3: Radial velocities of the WR and components of WR137. 
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Tabic 3: Radial velocities of the WR and components of WR137. 
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Tabic 3: Radial velocities of the WR and components of WR137. 
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Tabic 3: Radial velocities of the WR and components of WR137. 



HJD- RV HJD- RV HJD- RV 

2440000 km s' 1 2440000 km s' 1 2440000 km s' 1 

11389.68728 3lT~ 

11390.56354 -18.89 

11390.56831 -2.54 

11390.77701 -2.9 

11439.50374 -17.49 

11440.50850 -20.14 

11441.51545 -28.43 

11450.51896 -29.33 



This Table is the electronic version Table 3. n corresponds to the total 
number of spectra for one observatory. 
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